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Abstract 

Circadian clocks ubiquitous in life forms ranging bacteria to multi- 
cellular organisms, often exhibit intrinsic temperature compensation; the 
period of circadian oscillators is maintained constant over a range of phys- 
iological temperatures, despite the expected Arrhenius form for the reac- 
tion coefficient. Observations have shown that the amplitude of the oscil- 
lation depends on the temperature but the period does not — this suggests 
that although not every reaction step is temperature independent, the 
total system comprising several reactions still exhibits compensation. We 
present a general mechanism for such temperature compensation. Con- 
sider a system with multiple activation energy barriers for reactions, with 
a common enzyme shared across several reaction steps with a higher ac- 
tivation energy. These reaction steps rate-limit the cycle if the temper- 
ature is not high. If the total abundance of the enzyme is limited, the 
amount of free enzyme available to catalyze a specific reaction decreases 
as more substrates bind to common enzyme. We show that this change in 
free enzyme abundance compensate for the Arrhenius-type temperature 
dependence of the reaction coefficient. Taking the example of circadian 
clocks with cyanobacterial proteins KaiABC consisting of several phospho- 
rylation sites, we show that this temperature compensation mechanisms 
is indeed valid. Specifically, if the activation energy for phosphorylation 
is larger than that for dephosphorylation, competition for KaiA shared 
among the phosphorylation reactions leads to temperature compensation. 
Moreover, taking a simpler model, we demonstrate the generality of the 
proposed compensation mechanism, suggesting relevance not only to cir- 
cadian clocks but to other (bio)chemical oscillators as well. 



1 Introduction 

The circadian clock is one of the most remarkable cyclic behaviors ubiquitous 
to the known forms of life, ranging from the unicellular to the multicellular 
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level — including prokaryotes. Because of its importance, the underlying chemi- 
cal reactions have been the subject academic interest for a long time and have 
recently been elucidated experimentally. Circadian clocks have three important 
features: 

1. They persist in the absence of external cues with an approximately 24-h 
period, which is rather long compared with most chemical reactions. 

2. They can be reset by exposure to external stimuli such as changes in illumi- 
nation (dark/light) or temperature. 

3. The period of the circadian clock is robustly maintained across a range of 
physiological temperatures, (temperature compensation) pQ[2] . 

The emergence of cyclic behavior with a capacity for entrainment is theo- 
retically understood as being result of the existence of a limit-cycle attractor 
in a class of dynamical systems described by chemical rate equations. How- 
ever, the phenomenon of temperature compensation is not yet fully understood. 
Generally, the rate of chemical reactions depends strongly on the temperature. 
Most biochemical reactions, in particular, have an energy barrier that must be 
outcome with the aid of enzymes, and thus the rate can be expected to follow 
the Arrhenius form. Thus, the period of chemical or biochemical oscillators can 
be expected to strongly depend on the temperature f^- Thus, the ubiquitous 
temperature-compensation ability of biological circadian clocks, suggests that 
there may be some common mechanism(s) behind it. 

Overall, there are two possibilities: one is that compensation exists at each 
elementary step, and the other is that compensation occurs at the system level — 
for the total set of enzymatic reactions. Recently, it was found that the rates of 
some elementary reaction steps in circadian clocks depend only slightly on tem- 
perature [HIS] , suggesting that the activation energy barrier of some reactions is 
rather low. Although some element-level compensation is important, it is diffi- 
cult to imagine that every reaction step is fully temperature-compensated at the 
single-molecule level. Indeed, if that were the case, temperature changes would 
not influence the oscillation of chemical concentrations at all. However, even 
though the period of oscillation is generally insensitive to temperature changes, 
it is known that the amplitude changes indeed depend on temperature [5HH]- 
Furthermore, circadian clocks are known to entrain to external temperature cy- 
cles, and they can be reset by temperature cues [6l[9l[T0]. Thus, although tem- 
perature changes can influence the oscillation, the period is still robust. Hence, 
it is necessary to search for a general logic that underlies the temperature com- 
pensation phenomenon at the system level. Indeed, several models have been 
proposed [TT| - fT5] . In most of the studies, several processes that are responsible 
for the period are considered, which cancel out the temperature dependence 
with each other. With such balance mechanism, the temperature compensation 
is achieved. These mechanisms, however, need fine-tuned set of parameters, or 
rather ad-hoc combination of processes for the cancellation. Considering the 
ubiquity of temperature compensation, a generic and robust mechanism that 
dose not require tuning parameters is desirable. Here we propose such a mecha- 
nism that has general validity for any biochemical oscillator consisting of several 
reaction processes catalyzed by enzymes. The mechanism can be briefly outlined 
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as follows: 

Biochemical clocks comprise multiple processes, such as phosphorylation and 
dephosphorylation, with each generally having a different activation energy bar- 
rier IS.E. The rate of such reactions, then, is proportional to A/ exp(— Ai?/T), 
where Af as, the concentration of "free" enzyme available. If the enzyme con- 
centration were insensitive to temperature, the rate would just agree with the 
simple Arrhenius form, thus implying high temperature dependence. Now, con- 
sider a situation where several substrates share the same enzyme. At lower 
temperature, the reactions with higher activation energy will be slower and the 
substrates involved in these reactions will accumulate. Then, because they share 
the same enzyme, competition for the enzyme will also increase. Accordingly, 
the concentration of available (free) enzyme decreases, and when it reaches a 
level satisfying A/ exp(— Ai?/T) ^ 1, these enzymatic reactions will be highly 
suppressed. The system spends most of its time under such conditions, which 
limits the rate. Thus, the Arrhenius-type temperature dependence is compen- 
sated for by the concentration of available enzyme as Af exp{—AE/T) ~ 1. 

Although this is a rather basic description, its validity may suggest that 
temperature compensation emerges in any general chemical oscillation consist- 
ing of steps, catalyzed by a common enzyme. Here, we first study the validity 
of this enzyme-limited temperature compensation mechanism for the specific 
case of Kai protein clock model introduced by van Zon et al. [TB] to explain 
the circadian clock of KaiABC proteins in cyanobacteria, which was discovered 
by Kondo and his colleagues t8j[17| Indeed, in this system, the period of circa- 
dian oscillation of Kai proteins is temperature compensated |3- Further, some 
of the elementary components in this system, specifically, KaiC's ATPase ac- 
tivity and KaiC's phosphatase activity, were suggested to depend only slightly 
on temperature IHHB], but the origin of system-level temperature compensa- 
tion has not yet been explained. Here, we show numerically that system-level 
temperature compensation emerges from the differences in activation energy be- 
tween phosphorylation and dephosphorylation and competition for KaiA as the 
enzyme that catalyzes KaiC phosphorylation. Furthermore, we elucidate the 
conditions necessary for this temperature compensation to work. Then, based 
on this analysis, we introduce a simpler model consisting of few catalytic reac- 
tions, to illustrate the above mechanism. Possible relations between our results 
and reported experimental findings for circadian clocks are also discussed. 

2 Model 

The Kai-protein-based circadian clock, discovered by Kondo's group consists 
of KaiA, KaiB, KaiC proteins with ATP as an energy source [8]. KaiC has a 
hexameric structure with six monomers, each with two phosphorylation sites 
[19]. It has both self- kinase activity and self-phosphatase activity, but the self- 
phosphatase activity is usually stronger, and so it is spontaneously dephospho- 
rylated [201121] . KaiA, in a dimer form [19] , attaches to KaiC and thus increases 
its kinase activity, leading to phosphorylation of KaiC [2T1[22], while KaiB in- 
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hibits the activity of KaiA This phosphorylation/dephosphorylation 

process of the KaiABC proteins constitutes a circadian rhythm. Here, we sim- 
plify the process, to focus on the temperature compensation of the period. We 
reduce the two phosphorylation residues to just one, because abundance of singly 
phosphorylated KaiC is strongly correlated with that of doubly phosphorylated 
KaiC [23] so that the phosphorylation of the two residues are equilibrated on a 
rather short time scale. Next, we do not include KaiB explicitly in our model, 
because changes in the concentration of KaiB affect the period only slightly [33] . 
Note that although KaiB is necessary to generate circadian oscillations, the ef- 
fect can be accounted for by introducing a parameter value for KaiA activity. 
Here, we adopt a slightly simplified version of the model introduced by van Zon, 
et al. [16] (See also [26] ) 

First, each KaiC monomer has two states — active and inactive. Second, al- 
losterically regulated KaiC hexamers in the active state can be phosphorylated, 
whereas those in the inactive state can be dephosphorylated. A phosphorylated 
KaiC monomer energetically prefers the inactive state, whereas a dephospho- 
rylated KaiC has the opposite tendency. Here the flip-flop transition between 
active and inactive states occurs only from the fully phosphorylated or fully 
dephosphorylated states, as assumed in the concerted MWC model [17]. No 
intermediate states are assumed. Hence, the reaction process exhibits a cyclic 
structure as in Fig. 1. 



KaiA Active form 




Inactive form 



Fig. 1: Simplified model of circadian clock involving KaiAC proteins based on van Zon et al. [HI 
A KaiC hcxamer, with six phosphorylation sites, can also take an active or inactive state. In the 
active state, KaiC phosphorylation is catalyzed by KaiA at the rate fcp, while in the inactive state, 
dephosphorylation progresses without any enzymes at the rate k^p. Affinity between active KaiC 
and KaiA reduces as the number of phosphorylated sites of KaiC increases successively. 



Next, KaiA facilitates phosphorylation of active KaiC with an affinity that 
depends on the number of phosphorylated residues of each KaiC hexamer. 
KaiCs with a smaller phosphorylation number have stronger affinity to KaiA 
and are phosphorylated faster. This assumption is necessary for generating 
stable oscillations [TBI. Then, the reactions are given by 
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Ce A Ce (1) 

Co A Co (2) 

C. a-i (3) 

C, + A^=iAC, % C+i + A (4) 



Here Ci and C^ denote the concentrations of active KaiC and inactive KaiC, 
respectively, with i phosphorylated sites; A denotes the concentration of free 
KaiA dimer. To study the temperature compensation of the period, we must 
also account for the temperature dependence of the reaction rate. Here, the 
rates of phosphorylation and dephosphorylation are governed by the Arrhenius 
equation. Then the rate constants kdp and kp are as follows: 

kdp = Cdp e-Kp{-l3Edp) (5) 
kp = Cp exp{-PEp) (6) 



with inverse temperature (3 {(3 — 1/T) by taking the unit of Boltzmann con- 
stant as unity. We could include the temperature dependence of the rates / 
and & in a similar manner, but as the reaction between active and inactive 
states progresses faster and does not influence the period [TB], this dependence 
is neglected, van Zon et al. demonstrated that temperature compensation oc- 
curs when the speed of phosphorylation and dephosphorylation is completely 
temperature-compensated at a level of elementary reaction process. However, 
it cannot explain why the amplitude of oscillation depends on temperature, or 
the entrainment to temperature cycle occurs. The compensation mechanism at 
a system-level is wanted. Here, the formation and dissociation of KaiAC com- 
plexes occur at much faster rates than other reactions and so are eliminated 
adiabatically. Thus, the change in the concentration A is given by 

i=0 * 

where Ki{= kf^/k'^f) are the dissociation constants. Considering the increase 
in affinity for KaiA with the number of phosphorylated sites, we set Ki = 
Koa^{a > 1.0). We adopted the deterministic rate equation given by the 
mass-action kinetics, and it is simulated by using the fourth-order Runge-Kutta 
method. 
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3 Result 



3.1 Oscillation at varied temperature 

As mentioned, a KaiC allosteric model was analyzed. Specifically, we study the 
case in which the activation energy for phosphorylation, Ep, is larger than that 
of dephosphorylation, E^p- See the supplementary information for other cases. 
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Fig. 2: Oscillations in KaiC phosphorylation at various temperatures {0 — 1.0, 1.5, 2.0). (A) Red 
line indicates the time series of the mean phosphorylation level defined by ^^—QiCi/Ctotal , whereas 
the green line indicates that of the fraction of free KaiA, A/Atotal- A decrease in temperature 
causes a decrease in the amplitude of the phosphorylation level. (B) Time course of the abundance 
of each form of KaiC, Cj. At low temperature, the basal amount of C5 (magenta line) is remarkably 
high. 



For a certain range of parameter values, we found periodic oscillation in the 
KaiC phosphorylation level and free KaiA abundance, as shown in the time 
series in Fig. 2A. The oscillation is described by a limit-cycle attractor, as 
represented in the orbit in a two-dimensional plane of KaiC phosphorylation 
level and free KaiA abundance. As the temperature increases, the amplitude of 
the limit-cycle increases (Fig. 2A) . Lowering the temperature causes a decrease 
in the maximum amount of free KaiA and increase in the minimum level of 
KaiC phosphorylation. Further lowering it, however, results in the limit-cycle 
changing into a stable fixed point via Hopf bifurcation. The timeseries of Ci {i = 
0, 1, 2, .., 6) and the temperature dependence are shown in Fig. 2B for /3 — 1.0, 
1.5, 2.0. With a decrease in temperature, we see an increase in C5, that is. 
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the abundance of KaiC with five phosphorylated residues (5P-KaiC), which 
leads to an increase in the minimum KaiC phosophorylation level, as shown in 
Fig. 2B. Note that there is a remarkable change in the time course of C5 (the 
abundance of 5P-KaiC) at fi ~ 1.2 (Fig. 4A). Below this inverse temperature, 
the minimum C5 is close to zero. At higher /? (i.e., lower temperature), however, 
C5 never comes close to zero, and the minimum value increases with lowering 
of temperature. 

3.2 Transition to temperature-compensated phase 




Fig. 3: Dependence of amplitude and period of oscillation on the inverse temperature /? — 1/T. (A) 
Maximum and minimum values of mean phosphorylation level over a cycle. The maximum value is 
nearly constant against temperature changes, whereas the minimum value increases with ^ above 
/3c — 1/Tc ^ 1.2, i.e., at temperatures below the characteristic temperature Tc- The oscillation 
disappears via Hopf bifurcation at /3 J^i 3.5. (B) Period of the oscillation. At high temperature 
(low /3), the period changes exponentially with /3, whereas below Tc, the period is nearly constant 
against changes in temperature. 
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Fig. 4: Normalized concentration of each component, C^jCtotal (A), "^C/CtotaL (B), A/Atotal 
(C) plotted against the inverse temperature ^. (A)(B) Maximum (red), minimum (green), and 
average (blue) over a cycle. (A) Note that the maximum value of C5 is nearly constant, wheras its 
minimum and average increase with /3 beyond /3c (i.e., at lower temperature). (B) The average value 
of / Ctotal is fitted well by the value of the unstable fixed point of the equation (magenta line), 
which is proportional to cxp(— /3(i?p — E^p). (C) The concentration of free KaiA in the plateau 
region of C5 is plotted, as estimated from the time where C3 reaches a peak. This free KaiA 
concentration closely follows exp(/3i?p) (magenta line), for f3 > /3c, i.e., below the characteristic 
temperature. 



The transition at /3 ~ /3c = l/Tc ~ 1-2 is also reflected in the temperature 
dependence of the period. We plotted the period of oscillation as a function of 
(inverse) temperature, together with the maximum and minimum KaiC phos- 
phorylation levels (see Fig. 3). Above Tc, the temperature dependence of the 
period follows exp{(3Ep), as can be naturally expected from a reaction process 
with a jump beyond the energy barrier. However, at lower temperature, the 
period is no longer prolonged exponentially and is nearly constant. Thus, the 
temperature compensation of the circadian period appears at lower temperature. 
There is also a clear difference in the amplitude of oscillation below and above 
Tc- At higher temperature (without temperature compensation), the amplitude 
of oscillation is almost constant over a large interval of temperatures. However, 
at lower temperature (in the temperature-compensated phase), the amplitude 
decreases with lowering of temperature; eventually, the oscillation disappears 
via Hopf bifurcation. This decrease in amplitude is due to the increase in the 
minimum value of the KaiC phosphorylation level, caused by the increase in the 
minimum abundances of C5. The temperature dependences of the abundance 
of each KaiC and free KaiA also show distinct behaviors below and above 
(see Fig. 4, Supplementary Fig. 2). The average and minimum abundances 
of C5 increase remarkably with a decrease in temperature below Tc, whereas 
the maximum hardly changes with temperature (Fig. 4A). However, the am- 
plitude of the oscillation of Ci(0 < i < 4) has a peak at around Tc, and the 
minimum value increases as the temperature decreases. The minimum as well 
as the plateau value of free KaiA increase when the temperature is decreased. 
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The minimum inactive KaiC abimdance (independent of the residue number) 
and Cq showed behaviors similar to those of Ci(0 < i < 4), whereas the aver- 
age followed exp{—P{Ep — Edp)) dependence throughout the temperature range. 
Considering that the total amount of all KaiCs is conserved and the difference 
in activation energy, this dependence itself is rather natural. 

3.3 Proportion of free-energy of phosphorylation to that 
of dephosphorylation is critical for temperature com- 
pensation 

Thus, the transition in the oscillation and temperature compensation behavior 
at low temperature is the salient feature of the present system. Next, we ana- 
lyzed the conditions for Ep and Edp, the activation energies for phosphorylation 
and dephosphorylation, respectively, to elucidate the temperature compensa- 
tion. Supplementary Fig. 3 shows a plot for the region where temperature 
compensation appears in the parameter space of Ep and E^p. From Supple- 
mentary Fig. 3, we observe that the temperature compensation appears in 
the regime Ep > 5Edp- The overall periodic behavior is determined mainly 
by Edp/ Ep, rather than the individual magnitudes. When 0.2 < Edp/ Ep < 1, 
there still exists a transition to a phase with weaker temperature dependence of 
period on lowering the temperature, but this effect is not sufficient to produce 
the temperature compensation (see also Fig. 5A). 
(See supplementary informations for the case with Ep < Edp.) 
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3.4 Temperature compensation depends on KaiA and KaiC 
amounts 



A • Temperature compensated • No oscillation 
Phase transition • No transition 
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Fig. 5: Influence of activation energy difference and the KaiA abundance on the temperature de- 
pendence of the period. (A) Phase diagram of the temperature dependence of the cycle against the 
activation energy ratio {E^p/ Ep) and the abundance of total KaiA (A-totai)- Red: temperature- 
compensated oscillation with dr/d0 ~ (r: the period) below the characteristic temperature. 
Green: oscillation showing transition at the characteristic temperature. Blue: oscillation without 
transition, with a simple increase with temperature of the Arrhcnius form. Black: no oscillations at 
all. As \Ecip/ Ep\ increases from to 1, we successively sec temperature compensation, transition, 
and a disappearance of oscillation. (For the case with Ep < E^p, see supplementary information). 
As Atotai increases, the width of the temperature compensation region decreases. (B) Effect of 
KaiA increase on the period. An increase in KaiA led to narrowing of the range of the periodic 
solution, and also the range of temperatures at which temperature compensation occurs. 



As already mentioned by von Zon et al., oscillation of KaiC abundance requires 
that the amount of KaiA is less than that of KaiC [TB]. An increase in KaiA 
abundance leads to decrease in the period, finally leading to no oscillations. 
The range of temperatures where the oscillations exist narrows as KaiA abun- 
dance increases, and the oscillation disappears at higher temperatures. Here, 
the transition temperature Tc shifts to lower temperatures (see Fig. 5B). Fur- 
thermore, the temperature compensation at T < Tc is lost KaiA abundance 
is increased. Although the temperature dependence of the period is weaker at 
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T < Tc, the dependence exp{(3AE') still exists with AE' smaller than Ep. We 
plotted the range where temperature compensation occurs at T < Tc in the 
two-dimensional plane with KaiA abundances and Edp/Ep. We see that low 
level of KaiA abundance is necessary for temperature compensation. 

3.5 Range of temperatures where systems is temperature- 
compensated is narrowed by a reduction in phospho- 
rylation sites 

Since KaiC is a hexamer, we adopted six phosphorylation sites in our model. 
However, to understand the biological significance of this number of sites, we ex- 
amine models with a reduced number of phosphorylation sites (Supplementary 
Fig. 4) . We find that reducing the number of phosphorylation sites from hex- 
amer to pentamer, and then to tetramer, narrows the temperature range where 
oscillations exist and temperature compensation occurs. For the tetramer, tem- 
perature compensation is not observed even at T < Tc. 

3.6 KaiC phosphorylation cycle is entrained by tempera- 
ture cycle 



<j) = 71/2 n 

Temperature cycle 



3n/2 




Time (li) 



Fig. 6: Entrainment of oscillation to externally applied thermal cycles. Oscillations with various 
phases (0 — 0, 7r/2, tt, 37r/2) are entrained by the temperature cycle between /3 — 1.5 and 1.7 within 
about 10 cycles. 



It has been experimentally it is suggested that KaiC's phosphorylation cycle is 
entrained to external temperature cycle [lOj . To examine such entrainment, we 
cycle the temperature between /? = 1.5 and /3 = 1.7 (i.e., within the region in 
the temperature compensation) periodically in time with a period close to that 
of the Kai system (27 h). Within about 10 cycles, the phase of oscillation of 
the KaiC system is entrained with that of temperature, independently of the 
initial phase of oscillation. Thus, entrainment to temperature cycle is achieved 
(see Fig. 6). 
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4 Discussion 



Here we discuss how temperature compensation is achieved. As presented in 
the Results section, two stages are necessary: transition in the temperature 
dependence of the period, and complete temperature compensation of the period 
at lower temperature. As seen in the phase diagram (Fig. 5A), the former 
requires that Ep is sufficiently larger than E^p, which, as will be shown, means 
that the phosphorylation process is rate limited. For the latter, the abundance 
of KaiA should be sufficiently small so that there is limited free KaiA that can 
be used for phosphorylation (i.e., competition for free KaiA). As will be shown 
below, the abundance of free KaiA decreases as the temperature increases, which 
compensates for the increase in the rate constant of the reaction. 

1. Phase transition 

When there is a difference between the energy barriers for phosphorylation 
and dephosphorylation, the temperature dependence of the rate of each pro- 
cesses is different. Roughly speaking, the time scale for the phosphorylation 
process changes in proportion to Aexp{/3Ep), whereas the dephosphorylation 
has exp{(3Edp), where A is the concentration of free KaiA. Thus, there exists 
a characteristic temperature at which the two rates are comparable; the rate- 
limited reaction switches at this temperature, where the phase transition to 
temperature dependence occurs. Thus, 

^'-^ (8) 

where A, the concentration of free KaiA, is estimated from the steady-state 
solution of our model (see eq.d?])) to afford A ~ K^Atotai/K^ + C^ if Atotai <C 
and A ~ Atotai if Atotai > K^. 

Thus, a sufficient difference in activation energy is necessary for the phase 
transition-like behavior because the critical temperature will diverge as Ep c:^ 
Edp. If the difference is small, the critical temperature goes beyond the temper- 
ature for the onset of oscillation and the transition never occurs. Moreover, if 
Atotai is too large or too small, the critical temperature is lower or higher than 
the range where the oscillation exists, and thus for both the cases, the phase 
transition disappears. These estimates agree with the phase diagram in Fig. 
5A. 

2. Temperature compensation by self-adjustment of the KaiA concentration 
When the temperature is lower than the transition temperature, the phos- 
phorylation process takes more time, and C5 is accumulated before dephospho- 
rylation from Ce progresses, as already discussed (see Fig.2B). The increase in 
the abundance of active KaiC leads to competition for KaiA, and thus a decrease 
in free KaiA. If the total KaiA abundance is limited, the system reaches a stage 
where phosphorylation almost stops. This leads to the plateau in the time course 
of C5, as observed in Fig. 2B. This drastic slow down of the phosphorylation 
process occurs when A is decreased to the level at KpA < 1. Thus, during the 
plateau in C5, this approximate estimate gives, A oc l/Kp oc exp(/3Ep) , so that 
the temperature dependence of the phosphorylation rate KpA is compensated 



12 



for by the decrease in A. This plateau region is rate-hmited in the circadian 
cycle, making the whole period independent of temperature. 

In more detail, this compensation is also estimated as follows. The abun- 
dances of the inactive forms of KaiC decrease with an increase in C5. Consid- 
ering the differences in speed between phosphorylation and dephosphorylation, 
the total inactive KaiC abundance (at the fixed point) is estimated as 

SC* oc exp(-/3(^p - Edp) (9) 



which is consistent with Fig. 4. The flow from inactive KaiC is thus estimated 

by 

kdp^C* (xexp{-l3Ep). (10) 

This flow starts the phosphorylation processes from Co, but is slowed down 
at some residue number m. In the present model, this slow down starts at 
m « 2 ~ 3 due to the paucity of free KaiA. Following the above estimate of 
flow, the maximum Cm is estimated to be proportional to exp{—/3Ep). Now from 
eq.([71) and Km Ki {i > m) and Cj ^ Cm {j < m) at this time, the abundance 
of free KaiA can be estimated approximately as A ~ Atotai — (ACm)/ {Km + A) 
When Atotai is small enough, the minimum free KaiA abundance is smaller than 
Km,thnS A ~ Atotai - iACm)/{Km). 

A ^ , .r i-nr - — r ^MPEp). (11) 

Therefore, when Atotai is small enough and Ep is sufficiently greater than 
Edp^ the time scale of the phosphorylation process is kpA, and the period of the 
cycle is temperature-compensated. Indeed, as shown in Supplementary Fig. 2, 
the maximum Ci (and C2) values show exp(— /3ii^p) dependence, whereas free 
KaiA shows approximately exp(/3£'p) dependence when the phosphorylation is 
slowed down. 

In essence, the temperature compensation mechanism requires two proper- 
ties: 

difference in activation energy between phosphorylation and dephosphorylation 
processes, and a limited abundance of enzyme KaiA. The former is essential to 
the phase transition, and the latter to the compensation of the Arrhenius-type 
by the temperature dependence of the individual reactions. As long as these 
conditions are met, the period is temperature-compensated at low temperature, 
without the need for fine tuning of the parameters of the system. 

These two properties generally appear if there are two types of processes 
with different activation energy, with one type catalyzed by a common enzyme 
and the other without catalyzation. If the enzyme abundance is limited, the 
competition for the enzyme will lead to temperature compensation in cyclic 
reaction systems. In particular, if the main component of the chemical reac- 
tions has allosteric structure like KaiC, the competition for enzymes will occur 
naturally. 
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Fig. 7: Temperature compensation of period in a simpler model with enzyme competition. (A) 
Scheme of the simpler model. The substrate has five forms, and reactions So — > , Si — ^ 5*2 , 5'2 — > 
S3, S3 — ^ S4 are catalyzed by the same enzyme E, but the other reactions are not. The affinity 
between Si and E weaken as i increases. (B) Maximum and minimum values of S2 over a cycle 
and (C) the period, plotted against the inverse temperature (3. The oscillation is temperature- 
compensated over a specific range of temperatures, although the transition is not as sharp as in the 
Kai model. 



As an example, we consider the cyclic process shown in Fig. 7, where a cycUc 
change in five substrates occurs. Three processes — Sq Si ^ S2 ^ S3 — have 
a high activation energy barrier, and are catalyzed by the common enzyme E. 
In this case, the period of the oscillation in the concentration of each substrate 
is temperature-compensated at low temperature. (If the number of reactions 
with the common enzyme is decreased to two, compensation still appears, but it 
is weaker, as shown in Supplementary Fig. 5). This example demonstrates the 
generality of the present mechanism, and opens the possibility of applications 
to temperature compensation in other biochemical oscillations as well. 

The regulation of reaction process by autonomous changes in enzyme con- 
centration, as adopted here, was previously pointed out by Awazu and Kaneko, 
who reported that relaxation to equilibrium slowes down when the concentra- 
tions of substrate and enzymes are negatively correlated [28]. Excess substrate 
hinders the enzymatic reaction, leading to a plateau in relaxation dynamics. In 
the present model, the total concentration of KaiA as an enzyme is a conserved 
quantity but the fraction of free KaiA available is reduced when there is abun- 
dant active KaiC and so the speed of phosphorylation slows down dramatically. 
Such regulation of reaction speed by limitation of catalysts should be important, 
not only for the temperature compensation of circadian clocks in general, but 
also for other biological processes [29] . 

Now we briefly discuss the relevance of the present results to experiments 
on the Kai protein circadian system. In our model, it is assumed that KaiC's 
affinity to KaiA depends on the phosphorylation levels. This assumption is 
necessary not only for the emergence of temperature compensation but also for 
the existence of the oscillation itself. Recent reports that KaiC phosphorylation 
induces conformational changes may account for such changes in the binding 
affinity. 

In the actual circadian clock comprising Kai proteins, KaiB is also involved 
besides KaiA and KaiC. However, as mentioned, an increase in KaiB abun- 
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dance has minor influence on the clock [2S] . Earlier theoretical study suggested 
that KaiB binds to KaiAC complexes strongly and restricts the concentration 
of free KaiA [T^. Thus, the inclusion of KaiB is expected to not alter the 
present temperature compensation mechanism, but to facilitate it by strength- 
ening the limitation of enzyme availability. Note that previously reported tem- 
perature compensation at element-level — KaiC's ATPase activity [1] and auto- 
dephosphorylation activity [18^ is also relevant to the system-level compensation 
in our model. The former is relevant to achieve fast equilibration between KaiA 
and KaiA complex used for adiabatic elimination in eq.Q, and the latter to 
provide Edp < Ep. 

It is known that the amplitude of KaiC oscillation decreases as the temper- 
ature is lowered [8], which agrees well with our results (Fig. 2A). Indeed, as 
described already, this decrease in amplitude is tightly coupled with the tem- 
perature compensation mechanism. It is also interesting to note that in many 
circadian clocks, the period does not depend on the temperature, although the 
amplitude decreases with temperature |5]. If every elementary step of the cir- 
cadian clock were temperature compensated at the single-molecule level, this 
temperature dependence of the amplitude would not be possible. Furthermore, 
the entrainment of the KaiC oscillation to imposed temperature cycles, as ob- 
served in a recent experiment would also not be possible [TU] . In our study, the 
compensation is based on the different temperature dependences of the phos- 
phorylation and dcphosphorylation processes; the entrainment is a result of this 
difference. 

Our mechanism for temperature compensation depends on paucity of KaiA — 
increasing the KaiA concentration leads to loss of the compensation. We expect 
that this prediction will be directly confirmed in a future experiment. 
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Supplementary Fig 1: Temperature dependence of the amplitude and period of oscillation for Ep < 
Edp {Ep — 0.1, Edp — 1.0). (A) Maximum and minimum mean phosphorylation levels over a cycle. 
The maximum value is nearly constant against temperature changes. The minimum value increases 
with ^ for ^ > /3c — l/Tc 0.7, and is nearly constant below /3c- The oscillation disappears 
at 13 ^ 2.0. (B) Period of oscillation. The period changes with the Arrhcnius form, i.e., in an 
exponential manner over the whole range of temperature, while the slope changes below and above 
Tc. Below Tc (beyond 13c.), the temperature dependence of the period is about ex.p{l3Edp) , and is 
cxp{l3AE') with AE' larger than E^p at T > {/3 < pc)- 
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Supplementary Fig 2: Concentration of eaeli KaiC forms. Red line : maximum value, green line : 
minimum value, blue line : average value over a eycle. (A) Co, (B) Ci , (C) C2, (D) C3, (E) C4, 
(F) Ca, (G) Ce, (H) EC. 
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Supplementary Fig 3: Phase diagram showing temperature dependence of the period for various Ep 
and Efip, while Atotat is fixed at 1.2 x 10~^. Red: temperature-compensated oscillation satisfying 
dr/dji < (r : the period) below the characteristic temperature. Gray: oscillation without tem- 
perature compensation. When 0.2 < E^ip/Ep < 1, the oscillation is temperature-compensated at 
temperatures below the characteristic temperature. 
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Supplementary Fig 4: Temperature dependence of the period of reduced Kai models. The hexamcr 
is the original model shown in Fig. 1 containing Cq ~ Cq, while a pentamer and tctramer are the 
reduced models containing Co ~ C5 and Cq ~ C4, respectively. The same parameter values are 
used as in the original model. Models with fewer phosphorylation sites (z < 3) cannot generate 
any oscillation. In the reduced models, the temperature ranges where oscillations occur and where 
temperature-compensated oscillations occur are both narrower, and the temperature compensation 
in the tetramer model is not perfect. 
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Supplementary Fig 5: Temperature compensation in a simpler model with just two catalytic reac- 
tions sharing the same enzyme. (A) Scheme of the model. The reactions between Sq Si , Si 
S2-, S2 S3 arc catalyzed by the same enzyme E, and other reactions are not. The affinity between 
So and E is higher than that between 5*1 and E. (B) the period, plotted as a function of the inverse 
temperature 0. The range with temperature compensation is narrower, and the compensation is 
not perfect. (Inset: Maximum and minimum values of Si over a cycle). 
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